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Abstract 

We study the problem of calculating two-loop three-point diagrams with irre- 
ducible numerators (i.e. numerators which cannot be expressed in terms of the 
denominators). For the case of massless internal particles and arbitrary (off-shell) 
external momenta, exact results are obtained in terms of polylogarithms. We also 
consider the tensor decomposition of two-loop three-point diagrams, and show how 
it is connected with irreducible numerators. 
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1. The study of two- loop diagrams is required by increasing precision of experiments 
testing QCD and the Standard Model. As compared with the one-loop case, the calcula- 
tion of two-loop diagrams is technically much more complicated. Exact results are known 
for a few special cases only, mainly for two-point functions. There are, however, some 
approaches which make it possible to obtain numerical results [0, Q or analytic expansions 
providing reasonable precision for the cases of interest. 

Three-point two-loop diagrams are much less investigated. In this case, we have two 
independent external momenta (from which three independent invariants can be con- 
structed). Although some of the numerical approaches mentioned above (see also in |^) 
can also be applied to three-point functions, this usually requires a larger number of para- 
metric integrations to be done numerically, and the structure of the singularities in the 
integrand is more complicated. Asymptotic expansions can also be constructed, but they 
involve three- fold series in external momentum invariants (see, e.g., in |^). Furthermore, 
the expansion at large external momenta requires more information about three-point di- 
agrams with massless internal lines involving not only higher powers of the propagators, 
but also some numerators which cannot be cancelled against any of the denominators. We 
shall refer to these as "irreducible numerators" . One also obtains them in realistic calcu- 
lations with vector and spinor particles. Moreover, the problem of tensor decomposition 
of integrals with uncontracted Lorentz indices is also closely connected with irreducible 
numerators. This is one of the essential complications of two-loop three-point calculations 
as compared with both one- loop case p, and two-point self-energies 

The problem of irreducible numerators (and the related problem of tensor reduction) 
is the main subject to be discussed in the present paper. In fact, it can be considered 
as a continuation of the papers P, ^ where some exact results for three-point two-loop 
diagrams were derived^. The remainder of the paper is organized as follows. In section 2 
we collect some useful results for one-loop triangle diagrams. In section 3 we present some 
results for two-loop three-point functions and discuss the problem of irreducible numera- 
tors. Then, we calculate integrals with irreducible numerators for the planar (section 4) 
and non-planar (section 5) cases. In section 6 we consider the tensor decomposition of 
three-point two-loop diagrams. In section 7 we discuss the results. 

2. Here we shall briefly summarize some useful formulae for one- loop integrals. 
Definition. The n-dimensional one-loop three-point Feynman integral is defined as 



J{n; z/i,z/2,z/3| pI,pI,pI) = J 



where all external momenta Pi,P2,P3 (such that Pi + P2 + P3 = 0) are ingoing, and 
the "causal" prescription {q^)'" (g^ + iO)"'' is understood. Below we shall omit the 
arguments p\,p\,p\ in J . 

''Uniqueness" relations. At special values of the sum of the powers of denominators, 
X] z^j = z^i + z^2 + ^^3, the following formulae [0, [T^ are valid: 



J(n; 1^1, 1^2, ^^3; 



^n/2 -l-n 
n ' 1 



T{n/2 - u,) 1 



r(z/,) (p2)n/2-., 



^In the papers ml Ol some results were generalized to the case of an arbitrary number of loops. 
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I UiJ{n; z/i + 1, 1/2, z/g) + U2J{n; z/i, z/g + 1, z/g) + u-sJ^n; z/i, z/g, z/3 + 1) 



T,Ui=n—2 



n/2a-nA r(n/2-Z.,-l) 



Note that the sum of the powers of inverse momenta squared on the r.h.s. of is n/2. 
The case z^i = z/2 = z/3 = 1, n = 4. The result for this case is well-known (see, e.g., in 



15], |T^). Following the notation of refs. |]14[ g], we write it as 

CW{pIpIpI) ^ J(4; 1,1,1) = in' (p^)-i (4) 

x=pI/pI , y = pl/pl- (5) 

The function can be represented as a parametric integral, 

1 

^^'\x,y) = - f ^ —— (ln^ + 21ne), (6) 



(7) 



or in terms of dilogarithms, 

<|.(i)(x, = ^ I 2 (Li2 i-px) + Li2 i-py)) + In ^ In + In(px) \n{py) + y j 

-^(a;, = - a; - - 4:xy , p(x, = 2 (1 - x - ?/ + A)"\ (8) 

r/ie case z/i = z/2 = z/3 = 1, n = 4 — 25. In ref. 0, the following representation was 
obtained for arbitrary n = 4 — 2e : 

j(A-9r 1 1 7r^-i^+^- r(i + e)r^(i-e) 1 / d^rHiyQ-'-i^/Q-n 

^ '''^ iplV^' r(l-2.) ^/ (ye + (l-x-y)e + a;r^- 
The expansion in e yields 

J(4 - 2e- 1, 1, 1) = vr^- i^+2^(j9^)-i- r(l + e) {<!>^'\x, y) + £ ¥'\x, y) + ^(e^)} , (10) 

where is defined by (§p-([7p while can be represented as 

^) = - / + + . { ! + '"<^«' + <^-^-^'« + ^> 

-21n?/ln^-21n2^- -ln(a;2/)ln-|. (11) 

The result in terms of polylogarithms (up to the third order) is presented in eq. (29). 
The case z/i = z/2 = 1, z/3 = l-|-(5, n = 4. In ref. p[, the following formula was derived: 



J(4-l 1 l + ^) = ^l/df (^^)-' - (^/^)"' (12) 
The first terms of the expansion of J(4; 1, 1, 1 + 5) in 5, up to and including 5^, are 

1- ]-5{\Yix+\Yiy) + ^52Qn2 3,+ina,in^ + in2^)) $^^^2^,?/) + ]:6''^^'^\x,y)\{U) 
2 6 / 3 J 



iTT^ 



(Pi) 
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where is connected with the two-loop planar diagram (see in Section 3). 

Note that in all integral representations @, (H), (|TT]) and (|^) the denominator of the 
integrand can be represented as a propagator, p| (y^^ + (1 — x — y)C, + x) = {pi + ^^2)^- 



3 There are two basic "topologies" of the two-loop three-point diagrams: planar 
(Fig. la) and non-planar (Fig. lb) ones. Note that the diagram in Fig. lb is symmetric 
with respect to all three external lines, while the diagram in Fig. la is symmetric with 
respect to the two lower lines only. The corresponding Feynman integrals are: 



dq dr 



{pi + qY (pi + <? + ry {p2 - r)2 {p2 - q - r)^ g2 



(14) 



(15) 



Pi + q \ P2- q 



Pi + r 



P\ 




P2-r 



P2 



Pi + q + r \ p2 - q - r 



pi + q 



Pi 




P2- r 



P2 



(a) 



(b) 



Figure 1: Two-loop three-point diagrams: planar (a) and non-planar (b) cases 
The planar diagram (|lj) was calculated (in four dimensions) in ref. ^j. The result is 

c('\plplp^) = {i7rrip^)-'<^''\^,y), m 

where the function (see also (0)) is 



#P'(a:.a) 



It can be expressed in terms of polylogarithms up to the fourth order, see e.g. [|10 

eg. (19)B 

"^We would like to note that recently this result was checked numerically ipTf. 



As to the non-planar diagram ([T5|), in refs. ^S], ^ it was obtained (also at n = 4) that 



ITT 



'''' ipir h^'\x,y) 



with <l>(i)(x,|/) defined by (1)-®. Therefore, C^^) can be expressed in terms of diloga- 
rithms and their products. 

Now, let us consider the cases when some numerators occur in the integrands on r.h.s.'s 
of ([T5|) . First of all, let us introduce the following notation: 



C^'^^ [something] = {integral (14) with [something] in the numerator} , (19) 

C*-^-* [something] = {integral (15) with [something] in the numerator} . (20) 
In this notation, eqs. ( p!4D and (p^) correspond to C*-^-'[l] and C7*^^^[l], respectively. Inserting 



some numerators into (^If ) and (|l^) may produce divergent integrals. In these cases, 
we shall understand that dimensional regularization [^ is employed to regulate these 
singularities^, and (for two- loop integrals) we shall usually omit terms vanishing as e — >■ 0. 

In the paper [^ some integrals of the type of (p!9|), (pO|) were considered, and the results 
can be written (in the new notation) as 



pI c^'\pI,pIpI) 



ivr 



2\2 



Pi 



$(2) 



1 1 

x' X 



\{pi + q + rf] = \{P2 - q - rf] = pj C^'\pIpIpI) 



ivr 



2\2 



pI 



cm 
cm 



{q-rf 

{pi+qf 



i2+4e^4-2. YHl + e) 



(p2U+2e 



2e 



•2+4e^4-2e Y^il+e) 



{Pl 



2e 



(21) 

$(2)(a;,^), (22) 
(23) 



--^In^Uw C-A^^m ri,^)L(24) 



-$«(x,i/) + vI/«(x,i/)l, 



cm 



{Pi + rf{p2 - ry 



cm 



j2+4e^4-2£ 



{Pi + qf{.P2 - r) 

{ 1 



TT 



{pir 2(l-2e)(l-35) \e 



2 \\ix\ny + (I- X -y)^^^\x,y) - '—\ . (25) 



One should remember that "+(9(£:)" is understood in all r.h.s.'s. Some other results can 
be obtained from (|2T|)-([25|) by using symmetry properties, or they may correspond to 
propagator-type diagrams, for example (see pOU ) 



cm 



■m'Y ipi)-' 6C(3). 



(26) 



In all these examples (pl|)-(|26|) the numerators can be cancelled against the corre- 
sponding denominators, and (effectively) we obtain diagrams where some of the lines 



^For simplicity, we put the dimensional regularization scale parameter /io = 1. 
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are reduced to points. There exist, however, some scalar numerators which cannot be 
represented in terms of the denominators of the diagrams in Fig. la,b; namely: 



and 



C(') Uq + rf 



(27) 



Other irreducible cases can be related to these two ones. It is easy to see that, if we 
introduce auxiliary "forward scattering" four-point functions according to Fig. 2a and 
Fig. 2b, these numerators will correspond to the dashed lines (missing in Fig. la and 
Fig. lb). Calculation of these integrals (|27| ) will be considered in the next two sections. 




Figure 2: Auxiliary planar (a) and non-planar (b) four-point functions 



4. In this section we shall examine the planar diagram with irreducible numerator, 
C^^^ [g^]. It is easy to see that it can be written as 







(pi + q? 








\2 



- 2p^, [(g - 



(28) 



where only the first term on the r.h.s. is singular in e. All the terms on the r.h.s. are 
known (see eqs. (pSj), (pTf), (p6D), exept the last one, Pi C^"^^ [{q — r)^], which is convergent. 
Let us consider pi C^"^^ [(g — r)^] at n = 4 and introduce auxiliary analytic regular- 



ization by multiplying the integrand with [{pi + r)^] ^[{p2 



\2]5 



Since 

the integral is convergent, the result should correspond to the limit 5 — >■ 0. Then, let us 
consider g-integration and integrate by parts^ (see in pT|): 



Pi 



dg (g 



{Pi + qf {P2 - qy [(g - ry 



1-5 



26 



Pi 



dg 



d 



{pi + g)^ {P2 - qY dqf^ 



[{q 



\2-i5 



^The calculation of C^^) [q"^] can alternatively be considered in the coordinate space. 
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1 f dq{pi,pi + q) 1 f dq{pi,p2-q) , . 



6 J [{p, + (p, - q)2 [(g _ r^^s (p, + qf [{p^ - qYY [{q - vY] 

Now, let us study the contributions corresponding to the two integrals on the r.h.s. 
separately. Each of them is sungular in 5, but together they should give a finite result. 

In the first contribution, the sum of the powers of {p2 — t)^, {p2 — qY and (g — r)^ in 
the denominator is equal to two (half the space-time dimension if n = 4). So, one can try 
to use the "uniqueness" relation @) to transform the product of these propagators into a 
triangle. However, we need to introduce an additional regularization by multiplying the 
integrand with [{pi + r)^]~"^ [{p2 — qYY^ [{q — (otherwise we would obtain "bare" 

singularities in separate terms). Since the contribution we are considering is regular in 
5' (but not in 5), we need to remember that, after the calculation is performed, we have 
to let 5' — > first. Now, we perform several transformations by using the relations (0) 
and (|]) (and also trivial formulae for one-loop two-point functions), and we arrive at the 
following representation of this contribution: 

1 /■ /■ dg dr (pi,pi -|- q) 



6 J J [{pi + r)2]i+5+5' [{p, + g)2]2 [(p2 - r)2]i+<5 [{p2 - g)2]i+'5' [r2]i-5 [(g _ r] 



2]-5-S' 



^^Hplr-' ^-(l + 5>? J(4;1,2 + 5M 



26^(1 + 6) I '"^^^ 6 



+j;iplY^' {(S + 5')il + S + S') J(4; 1,2 + 5 + 5',!)- 5(1 + 5) J(4; 1,2 + 5, 1)] j , (30) 

where the integrals J are one-loop triangles defined by (0). 

Considering the second term on the r.h.s. of eq. (|29| ) and using the fact that, due to 
momentum conservation, —{pi,P2 — l) = {P2,P2 — q) + {P3,P2 ~ q): it is easy to see that 
the term {p2,P2 — q) should give the same as ( pOD with pi p2. To calculate the term 
corresponding to {j>z,P2 — q), we need to use only the relation (H) (and we do not need to 
introduce additional parameter 5'). In such a way, we get 

Iff _ dg dr {p3,p2 - q) 



5 J J [{pi + r)2]i+5 (pi + g)2 [(p2 - r)2]i+5 [{p2 - g)2]2 [r^f-i [(g - rf] 
{J(4;l,l,l) + (p2)^ J(4;l + 5,1,1)} 



25(1 + 5) 

Iff dg dr 



2 ^ ^ [{pi + g)2]i+^ [{j>2 - r)2]i+'5 [{j,2 - g)2]i-5 [r2]i-5 (g - rf 



(31) 



Note that in the last integral on the r.h.s. the denominator (pi + r)^ is missing. As 5 ^ 0, 
this integral gives C*-^-* [{pi + r)^] which cancels the corresponding term in (p^). 

Careful analysis of the contributions (|30D , (|3l|) shows that what we need, in addition 
to the formulae presented in section 2, is the expansion of the integral J(4; 1, 1, 2 + 5) up 
to (and including) 5^ terms. To get it, it is convenient to use the representation (|^) with 
5 substituted by 1 + 5. The result of this calculation can be presented as 



J(4; 1,1,2 + *) = J^,^ TT^ {"I ''""^ + I ^ + In^l'iS + 1" 
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H — (In a; + In x In ?/ + In a; In y + In y 
6 ^ 

[\ - ^(Inx + Iny)^ (1 - X - ^^^\x,y) + ^ fi^^H^^,?/) + OiS"^ .(32) 

Here we define a set of functions fi^^^ via the derivatives of the functions (the general 
case of these functions was considered in [|1^]) as 



fiW(x,y) = A 



5a; ^ ' dy 

For = 1, the corresponding function fi^^-* is trivial, 

y 

Vl'^^'' (x, y) = In a; + In y — (x — In 



X 



For N = 2, the function is defined by eq. (|17|) , and fi'-^'' can be represented as 



n^^\x,y) = ^Inxln^^ + S / ^ (ln^ + 21ne)ln 





(1 - X-?/)^ + x' 



X 



(33) 



(34) 



(35) 



or in terms of polylogarithms, 

(a;, y) = 6 [Li3 i-px) + Li3 (-P2/)] + 3 In ^ [Li^ {-px) - Lis (-P2/)] 

X 

~ In^ ^ [ln(l + px) + ln(l + py)] + ^(71^ + In(px) ln(/«/)) [In(px) + \n{py)] . (36) 
Note that, combining the derivatives of ^^^^ in a different way, we may obtain ^'•^^ as 



In^ ^^'\x,y) = - 
X A 



d d 
[l-x+y)x—[X^^^\x,y)) - {l+x-y)y— [X^'-^^x^y) 



dx ^ ' ' dy 

Finally, using (|3^) , the integral with irreducible numerator can be represented as 



(37) 



C(2) 



2\2 



A pi 

2\2 



[ITT 



2 Pi 



In X + In y — (x — y) In 

X 1 



X. 

y 1 



<l>«(x,y) 
fO(£). (38) 



Note that the factor multiplying is proportional to ^'•^•'(x, y), eq. 



5. Let us consider the non-planar diagram with irreducible numerator, C*^^-* [(g + r)^]. 
It is easy to see that it is convergent and can be written as 

^(2) [(g + rf]= pIC^^) [1] + C'(^) [(p2 - g - r)=^] - 2^(2) [(p^, - g - r)] . (39) 
Now, using the symmetry of the non-planar diagram it is possible to see thatQ 

C^^^ [{p2.P2- q -r)]= ^(2) [(p2, pi + g + r)] . (40) 

''This can be shown using the fact that the non-planar integral ( [1^ ) remains invariant when we simul- 
taneously change r ^ p2 — r and q —pi — q. 
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Writing the l.h.s. of ( ^0]) as half the sum of the l.h.s. and the r.h.s., we get 



'2,P2 -q-r)] 



Up2,Pi+P2) C^'^ [1] 



and the resuh for the considered integral reduces to 



- q- r)^] - (piP2)C'^^^ [1] • 



(41) 



(42) 



According to (|T]) (see also [|), C(2) [{p^ — g — r)2] can be represented in terms of the 
planar diagram as C^^-* [1]. So, finally we arrive at the following result: 



\{q + rf] =pI c^'\pIpIpI) - (pm) c^'\pIpIpI) 



(43) 



where the integrals contributing to the r.h.s. are defined by eqs. ( p!4D and (p^ . This 
equation is illustrated by Fig. 3 where "—1" means that the irreducible numerator can be 
considered as a negative power of the corresponding denominator. It is interesting that 
we have a "mixture" of different topologies on the r.h.s. of eq. (H3|) and Fig. 3. 




Figure 3: Result for the non-planar diagram with irreducible numerator 



6. In this section we shall demonstrate that the problem of tensor decomposition of 
two-loop three-point functions also requires integrals with irreducible numerators. 
Let us consider the case with one free Lorentz index. It is easy to show [^2|, Q that 

[q,] = A-'{{pIp,^ - {P^P2)P2,) [ip,q)] + {pIp2, - iPiP2)pi,) C^'^ [(^2?)]} (44) 



where A is nothing but the Kallen function, A 



pIpI 



' {P1P2) 



{j>1Y\^{x,y) . The same 



as (1^) is valid for C*-^-' [g^] (we only need to change C*^^^ — ^ (5*^^^ on the r.h.s.). Moreover, 
the results for C*^^-' [r^] and (5*^^^ [r^] can be obtained by changing q into r everywhere. 

For the two-loop planar diagram, we express {piq) and {p2q) in terms of the denomi- 
nators of the diagram in Fig. 2a, and we get irreducible numerators (|27|) , 



C^'^ [{piq)] = \ [-pIC^''^ [1] + C^'^ [{Pi + qf] - C^'^ [q'] } 



(45) 
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(46) 



Note that in C-^^ [(ps?)] the irreducible integrals C^"^^ [q^] disappear. In the case of C*^^^ [r^], 
we do not have irreducible numerators either, because q^ in (|45|)-(^6D should be replaced 
by r^, and the corresponding integral (^) is trivial. 

For the non-planar case, we need to express {piq) and {p2q) in terms of denominators 
of the crossed diagram shown in Fig. 2b. So, we get 



C^'^ [{Piq)] = I {-pIC^'^ [1] + [iPi + qf] - [i] } 



(47) 



{p2-q-rf 




{p2-rY 




\2 


+ ^(2) 


{q+rf 



and we see that the irreducible numerator (g + r)^ appears in (0). Analogous results for 
(7(2) [(p^r)] and C*^^) [(^2^^)] can be obtained from (^)-(|48|) by using the symmetry of the 
non-planar diagram. 

When we consider more free Lorentz indices, C^^^ [q^iqa] ,0^"^^ [q^r^], etc., we need in- 
tegrals with higher powers of irreducible numerators. Such integrals will be considered in 



more detail in 23 



7. In the present paper we have considered an approach to the calculation of three- 
point two-loop diagrams (Fig. la,b) with irreducible numerators (|27| ) which correspond 
to the dashed lines of auxiliary four-point diagrams shown in Fig. 2a,b. We have obtained 
exact results for both planar (Fig. la) and non-planar (Fig. lb) cases. The result for the 
non-planar case can be expressed in terms of the functions corresponding to the diagrams 
without numerators, as shown in Fig. 3. At the same time, the planar case is more 
complicated, and the result (|38D contains a new function n^'^^ which can be represented 
in terms of the derivatives of the function (see (PB|)). So, the full set of "non-trivial" 
functions occurring in two-loop three-point calculations (with massless internal particles 
and off-shell external momenta) are: 

eqs. (|^)-(|3)) corresponding to the one-loop triangle in four dimensions (H); 

eq. (|TT]) (and eq. (29) of |^), e-part of the one- loop triangle; 

eqs. (p^) (and eq. (19) of [l^), corresponding to the two-loop planar diagram; 

eqs. (|35|) -(p6D, arising in the two- loop planar diagram with irreducible numerator. 
It is interesting that all these functions can be obtained from expansions of dimensionally 
or analytically regularized one-loop integrals in the regulating parameters {e or 6), see 
eqs. ([T0|), (p!3|), (|32[). All of them can be expressed in terms of polylogarithms. Note that 
the results for some other diagrams can also be expressed in terms of these functions, for 
example: one- and two-loop off-shell four-point functions with massless internal lines [Q, 
one-loop diagrams with higher number of external lines [^, two- loop massive vacuum 
diagrams [^, ^ . 



We have also discussed the problem of tensor reduction of two-loop three-point inte- 
grals which also requires integrals with irreducible numerators. A recursive procedure for 
calculating integrals with higher powers of irreducible numerators can be constructed ||23|| . 
Such integrals are needed for the decomposition of tensor integrals with a larger number 
of uncontracted Lorentz indices. 
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Important applications of the obtained results are connected with the large momentum 
expansion of two-loop three-point massive diagrams, and with the calculation of vertex 
corrections in massless QCD. 
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